AF Inv. No. AFB00657 
High Resolution SAR Processing Using Stepped-Frequency Chirp Waveform 

STATEMENT OF GOVERNMENT INTEREST 

The invention described herein may be manufactured and used by or for the 
Government for govermnental purposes without the payment of any royalty thereon. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to radar systems and more specifically to a 
waveform process for use with radar systems. 

A Synthetic Aperture Radar (SAR) system is used for ground mapping even in 
poor weather conditions in which optical sensor cannot be used. Mapping range of SAR 
is generally longer than optical sensor and it can also be operated during night. Using 
SAR, desired resolution in azimuth is achieved by maintaining phase coherence between 
pulses and combining them in proper manner. The length of collected data (array length) 
is determined by desired azimuth resolution and collection geometry. 

The following references disclose recent SAR systems that can be used Avith the 
preset invention and the disclosures of which are incorporated ho-ein by reference: 

U.S. Patent No. 6,492,932, SYSTEM AND METHOD FOR PROCESSING 
SQUINT MAPPED SYNTHETIC APERTURE RADAR DATA, 
invented by Jin, Michael 

U.S. Patent No. 6,255,981, METHOD FOR RANGE ALLIGNMENT AND 
ROTATION CORRECTION OF A HIGH RESOLUTION IMAGE IN 
AN INVERSE SYNTHETIC APERTURE RADAR SYSTEM, invented 
by Samaniego, Raymond 
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U.S. Patent No. 5,430,445, SYNTHETIC APERTURE RADAR GUmANCE 
SYSTEM AND METHOD OF OPERATING SAME, Peregrim and 

U.S. Patent No. 5,424,742, SYNTHETIC APERTURE RADAR GUIDANCE 
SYSTEM AND METHOD OF OPERATING SAME, Long. Albert H. 

The above-cited references all describe SAR systems developed by Raytheon and 
which can use the present invention. The best reference is U.S. Patent No. 5,430,445, 
SYNTHETIC APERTURE RADAR GUIDANCE SYSTEM AND METHOD OF 
OPERATING SAME, Peregrim, and this system uses a chirped signal provided by a 
chirp generator, and its performance can be improved, as discussed below. 
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SUMMARY OF THE INVENTION 

The present invention is a waveform process for use with S AR groundmapping 
radar systems, and is based upon the following observations. 

Range resolution in SAR image is inversely proportional to the transmitted 
signal bandwidth in nominal SAR systems. Since there is a limit in the transmitted 
bandwidth that can be supported by the radar hardware, there is a limit in range resolution 
that can be achieved by processing SAR data in conventional manner. However, if the 
frequency band of the transmitted signal is skipped within a group of sub-pulses and 
received signal is properly combined, the composite signal has effectively increased the 
bandwidth and hence improvement in range resolution can be achieved. The proposed 
new and practical approach can effectively extend the lunit in range resolution beyond 
the level that is set by the radar hardware units when conventional method is used. 

Several people proposed some step-frequency schemes for high resolution SAR in 
the past. They all tried to combine data before dechirp processing and it is not considered 
to be practical for general high resolution SAR systems since dechirp ftmction for high 
resolution mode is performed before A/D conversion. The new method proposes range 
samples from sub-pulses to be combined after dechirp and A/D fimction and therefore it 
is easy to implement. 

Based on a clear understanding of signal relationship between sub-band steps it 
was necessary to derive the condition to connect returned signals from sub-pulses without 
discontinuity around the sub-pulse boimdaries. To maintain the phase coherence it was 
found that the proper phase adjustment and time shift were required. The idea was proved 
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valid through test by implementing it to the SAR simulation program. Final image and 
intermediate results conformed to the prediction. 

This invention claims that the foUowings are new. 

1 . The idea of combining sub-pulse range samples after the dechiip fimction 

2. Derived requirements of phase adjustment and time shift to maintain 
continuity between sub-pulse data 

3. Two alternative embodiments based on fimctional sequence for data 
composition and any version that are alternative but equivalent 
implementations 

4. Data cutting and concatenating to minimize the error effect around the step 
boundaries 

5. SAR groundmapping by transmitting stepped frequency chirp pulse 
waveforms with pulse pairs where the second chirped pulse has a higher 
frequency than the first chirped piilse. 
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DESCRIPTION OF THE DRAWINGS 



Figure 1 is a chart of transmitted pulses with instantaneous frequencies of two 
steps for increased bandwidth. Sub-pulse pair: (0, 1), (2,3), (4,5), ... ; 

Figure 2 is a chart of the frequencies of the returned signals in 3 steps; 

Figure 3 is a chart of signals from two subpulses with proper phase adjustments 
and time shift; 

Figure 4 is a block diagram of a range sample combining block in the SAR signal 
processing chain; 

Figure 5 is a chart of the data combining sequence of configuration A; 
Figure 6 is a chart of the data combining sequence of configuration B; 
Figure 7a-7d are charts of the formed image and range IPR plots with reduced 

single step bandwidth (a,c) and fiiU composite bandwidth (b,d); and 
Figure 8 is a chart of composite range samples by cutting arid concatenating 

dechirped range signals fiiom two sub-pvdse returns 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 



The present invention includes process for use with SAR, systems that uses a stepped- 
frequency chirp waveform during groundmapping. 

Range resolution in SAR images is determined by the radar bandwidth: range resolution 
is inversely proportional to the radar bandwidth. Therefore, for very high resolution in range, 
radar bandwidth of the SAR system is required to be very large. To meet the requirement of 
desired minimum average power with limited peak power, stretched pulse with linear frequency 
modulation (LFM) is widely used in SAR systems. As the radar bandwidth increases, hardware 
performance degrades and there is a limit on bandwidth that radar hardware such as antenna and 
receiver can support. Natural approach under this situation where very high range resolution 
requirement cannot be met due to hardware limitation is to transmit sub-pulses that have 
varying center frequencies with limited bandwidth and combine the return for each sub-pulse so 
that increased bandwidth is effectively achieved. This technique of transmitting sub-pulses with 
varying center frequencies and combining returned signals from sub-pulses to achieve 
effectively increased bandwidth is called stepped frequency or step stretch method. Fig. 1 
shows the transmitted pulses with two frequency steps to achieve wider bandwidth by 
combining returned signal from two sub-pulses. 

There were ideas proposed in the past on stepped frequency technique by combined chirping 
signal. Unlike the past work, the proposed method shows how to combine returned data from 
each sub-pulse after dechirp processing. It is considered to be practical and useful sine most 
high-resolution SAR systems perform range dechirp 
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processing before analog-to-digital (A/D) conversion and dechirped digital data is 
available at a rate that can readily be handled. Following discussions will provide detailed 
expressions to combine sub-pulse returns and other relevant techniques. 

To combine returned signals from sub-pulses with different center frequencies, it 
is important to imderstand the relationship between dechirped signal from each sub-pulse. 
Comparing the dechirped signal expressions from each pulse will lead to the adjustments 
that are required when sub-pulse returns are combined. Chirp slope y and pulse index n 
will be used in the signal expressions. Also, for generalized derivation with multiple 
frequency steps, step index m within each sub-pulse group will be used. Time delay to a 
target and reference point at step m will be denoted by tum and ts^m , respectively. Phase 
for transmitted and received signals and reference signal for LFM can be expressed as 
follows. 

Transmit: {t, n)^fmt-^^{t- tnf 

Receive: <Df (/,«) = fm{t-Tt,m) + ^{t-tn-Tt.mf 

Reference: Of'^C^/i) = /^^+^(/-f„-T.,m)'' 

Dechirp processing mixes the received signal with the reference signal to provide phase 
expression as 

<S>m{t, n) = + = fmTs. « + [/m + rO - Ts, m)](r/. « - r.. «) — ^ (Tt, m " Ts. mf (1) 

where new time variable i -t-tn. 

Let's consider a sub-pulse group that includes pulses with index m and m4. It is 
generally true that 
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However, it is important to note that difference between the time delay to a target and the 
time delay to the reference point stays about the same in two sub-pulses. That is: 

With this good approximation on time delay difference within the same sub-pulse group, 
we can rewrite (1) for reference step M and any other step m within the same sub-pulse 
group with the substitution of Ar = r/, m - r*. m as follows. 

<Dm ) = fmTs. «i [^w + /(^ - m)] A - At^^ 

To get a composite signal without discontinuity from sub-pulse retums it is required to 
maintain phase coherence. This means that the phase expressed in the above expressions 
should be equal between step M and m for all values of time variables i . By comparing 
phase terms that are independent of targets, required phase correction to the signal at step 
m is obtained as 

AOm = fidts^M — fmXs^m (2) 

To obtain the time delay adjustment at step /w, we introduce required time shifting 
variable Arm and set the multiplying terms to Ar/ in two expressions to be equal as 

^ + ;'(f-ATm-rj.m) 
Arm is then obtained as 

/lirm = Ts.M-Ts,m-{fM- fm)l Y (3) 
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Using the above expressions, phase record at step m can be adjusted for phase 
coherence to the reference phase record at step M through time shifting and phase 
adjusting as follows. 

Derived time shifting that is required to maintain phase continuity at step m can 
be validated by observing the relationship between neighboring steps illustrated in Fig. 2. 
Depicted in the figure are the frequencies of transmitted signals at three steps and their 
returns. Two groups of dotted lines indicate returned signals at two sub-pulses around the 
central one that is assumed to be the reference step. To maintain phase continuity, 
adjustments should be performed so that linear frequency profiles from all steps are 
connected. To find the required adjustments at step M-l to make it aligned to step Af, we 
use three points A, B, and C as marked in the figure. A and 5, respectively, are the time 
delays to the reference point before and after adjustments at step MA, and C is the time 
delay to the reference point at step M. It can be noted that time delays for points B and C 
are a/ - 1 and Ts, m , respectively. Then, the required time shift is BA and it cm be 
obtained as follows. 

Atm-} = BA = CA-'CB 

. ' This expression is identical to -(3)' with m-M-L 

The paragraphs above described requirements on phase and time shifts 
adjustments when received signals from sub-pulses that have different center frequencies 
need to be combined for the improvement of range resolution. In implementing the 
derived condition it is important to determine how we handle the data around the edges of 
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each band. Because of the inaccuracies around band edges and in the transition band it ii 
difficult to get good results by blending neighborhood sub-pulse data. A more desirable 
approach is cutting out data at the mid-point between neighboring sub-pulse returns and 
concatenating them. Another issue is to determine at what stage in the processing chain 
the sub-pulse returns should be combined. 

Let's consider the returned signals from two sub-pulses after proper phase and 
time shift adjustments as explained in the previous section. Fig. 3 illustrates time history 
of signals returned from near and far range edges of a swath through dechirp and range 
deskew processing that are normally required for high-resolution S AR systems with 
single-step frequency band. As observed in the figure, length of the composite signal 
from sub-pulses is increased and this leads to the improvement in range resolution. Since 
sub-pulse data are combined, the resulting PRF is reduced by a factor that is equal to the 
number of sub-pulses. The final effective PRF has a lower limit that needs to avoid the 
Doppler ambiguity. The lower limit requirement on PRF to avoid the Doppler ambiguity 
sets the upper limit for the time interval between sub-pulses. If the time interval is too 
small, target signal can significantly be corrupted by the target returns at different ranges 
from other neighboring sub-pulses. This effect due to increased PRF that is called the 
range ambiguity should also be avoided. To limit the range ambiguity effect it is 
desirable to maintain a small as possible overlap between each sub-pulse step. However, 
considering the uncertainties on band edges and effect of transition band, spectrum 
overlap should be allowed up to a certain level. Also, in some cases as discussed for the 
second configuration that will be presented, some overlap is required. 
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The next question is how the data in the overlapped portion are combined. 
Without accurate characterization of pass band and transition, blending data in the 
overlapped portion leads to an unacceptable level of error around the boundary between 
sub-pulses. It is believed that the best way of combining sub-pulse data is cutting the data 
in the middle of the overlapped region and connecting them since slight sample shift has 
less impact than phase or magnitude error on the performance. This cutting and 
concatenating technique will be used for two configurations described in this section. 

Required bandwidth and center frequency for each step will be calculated. Denote 
the band overlap normalized by single step bandwidth by . Then, the required fuU 
bandwidth BWwA single step bandwidth BW^ can be calculated as follows for the 
desired resolution 5 , number of steps AT ^ , mainlobe-broadening factor k, overlap , 
and speed of light C. 

kC 



BW = - 



2S 



BW,= (4) 

Center fi^quency at step m, for in=l,2,. . ., N^^ , is calculated by 



BW 

Um) = f,-^ + BW, 



(7«-l)(l-OJ + l 



(5) 



for the center fi-equency of the composite waveform, 

Two configurations proposed in the following for the implementation of 
combining dechirped return signals fcom multiple sub-pulses correspond to the common 
block marked by 'Step Stretch Configuration' in Fig. 4. 
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as 

since 



Return signals from sub-pulses can be combined after range-deskew processing 
shown in Fig, 5. The advantage of this approach is minimized spectrum overlap 
range samples returned from all sub-pulses are aligned. As mentioned, reduced spectrum 
overlap results in reduced range ambiguity. A drawback of this method is increased 
ripples around the sub-pulse boundaries caused by the range deskew function. Ripples 
may be able to be compensated up to a certain level using the predicted ripple output that 
can be calculated with ideal signal pattern. Discrepancy between the real and the actual 
signal shapes can result in ripples of unacceptable level. 

Required time shift to maintain continuity was derived above. Since the calculated 
time shift is not an integer in general, nearest integer value can be used to align range 
samples. Fractional shift can be implemented by applying linear phase after taking FFT. 
However, realization of fractional time shift this way can cause signal distortion as 
confirmed by simulation and therefore, is not recommended. Since absence of fractional 
timing shift has insignificant effect compared to the phase or magnitude error, in 
particular for targets near map center, just taking the nearest integer for the required time 
shift will lead to acceptable performance. 

An alternative approach for combining range samples from multiple sub-pulses is 
to perform data combining fimction before range deskew as shown in Fig. 6. With this 
fimctional sequence, there are no large ripples that occur around the sub-pulse boundaries 
with the configuration A since range deskew is performed after the sub-pulse returns 
combined. A drawback of this approach is the increased overiap between sub-pulses to 
insure no data gap caused by cutting data in each pulse return. Since time delay 
corresponding to swath width is generally small compared with the pulse width in most 
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high resolution SAR systems, required overlap bandwidth is generally small compared 
with the radar bandwidth. As observed from tfie data shape in Fig. 6, required minimum 
overlap between sub-pulses is determined by swath width. For swath width SWand pulse 
width minimum overlap normalized by the pulse width is calculated as 
2^ 

■ cr. 

Overlap bandwidth is then computed by 

over C T 

where BW^ is the bandwidth of the single step expressed in (4). 



SIMULATION RESULTS 

Test target signal could be generated and processed using a simulation program. 
Using a single step with reduced bandwidth, formed image shows degraded range 
resolution in Fig. 7(a). Combining two sub-pulse retums using two-step radar frequency 
bands, range resolution of target image is improved as shown in Fig. 7(b). Impulse 
response in range in each case is plotted in Fig. 7(c) and 7(d), respectively. Since the 
target image is displaced in range from map center, dechiiped signal should be sinusoidal 
with non-zero frequency. The real part of the composite signal based on configuration B 
is illustrated in Fig. 8. Plot of the composite signal shows smooth connection between 
retums from two sub-pulses. 

The present invention uses the terrain classification systems previously developed 
by Raytheon, as described in the Peregrim patent, for ground mapping. As discussed 
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therein, the feature classification process includes, in more detail, the operator classifying 
the features according to the typical feature types described in Table 1. 



Table 1. 



NATURAL 


MAN-MADE STRUCTURES 


Concrete 


Simple Roof Building 


Asphalt 


Complex Face Structure 


Grass 


Complex Roof Structure 


Trees 


Complex Structure 


Dirt/Gravel 


Bridge 


Water 


Complex Volunie 



The latter is done by selecting a feature type from Table 1 . and placing the cursor 
inside each of the homogenous regions thus labeling each pixel in the homogenous region 
as such as feature. Natural features are classified based on texture, gray level and context. 
Man-made structures are identified from shadows and line of sight geometry and are 
classified based on texture and context. Man-made structures are classified as follows. 

If a top of a structure has a feature which is smooth, it is classified as a simple 
roof building. Typically, the leading edge dihedral established by the vertical surface of 
the building and the ground provides the radar reflectivity. The radar cross section (RCS) 
of the top of the building, the trailing edge of the building and the ground is lower than 
the leading edge dihedral. Aircraft hangers and POL tanks with smooth tops also exhibit 
such radar reflectivity. 
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A complex face stracture has a feature that has reflectivity along the entire side of 
the structure. An example of this type is a building having windows along the entire 
vertical face. In order to identify and classify a building as a complex face structure, it is 
necessary to have oblique photography. 

If the top of a structure is rough in texture due to, for example, air conditioning 
and plumbing equipment fixtures, it is classified as a complex roof building. The entire 
roof will exhibit a high RCS due to scattering effects fi-om the fixtures. Power plants 
exhibit such radar reflectivity. 

A feature is a complex structure when the feature exhibits reflectivity over the 
entire surface of the structure. This can be represented by a combination of a complex 
roof and a complex face structure. Refineries and metal pipelines typically have features 
with these characteristics. 

A feature is classified as a bridge if it is a slab bridge or a road bridge. In general, 
this category includes all bridge types except for truss bridges which exhibit a higher 
radar cross section. 

Complex volume is a structure with features that exhibit reflectivity over the 
entire volume of the structure. Examples of this type of structure include truss bridges, 
transformer yards and substations. 

Each pixel in each like region is coded to identify the feature type (i.e. here codes 
are used which corresponds to the features of Table 1 .). At a boundary, the code of each 
pixel on each side of the boundary indicates the type of feature pair (i.e. a code pair 
indicates a grass/water feature pair). 
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Once the operator has classified the features of the structures within the template, 
these regions are specially coded with respect to their class types. 
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